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In this paper, we use the strong held limit approach to investigate the gravitational lensing 
properties of braneworld black holes. Applying this method to the supermassive black hole at the 
centre of our galaxy, the lensing observables for some candidate braneworld black hole metrics are 
compared with those for the standard Schwarzschild case. It is found that braneworld black holes 
could have significantly different observational signatures to the Schwarzschild black hole. 
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I. INTRODUCTION 

The braneworld paradigm provides an interesting 
framework within which to explore the possibility that 
our Universe lives in a fundamentally higher dimensional 
spacetime. Unlike the Kaluza-Klein picture where the 
extra dimensions must be compactified on a length scale 
R < 10~ 17 cm in order to evade our detection, confine- 
ment of the standard model fields to a 3-brane, with only 
gravity propagating in the bulk, allows large, and even 
infinite, extra dimensions. 

Recent work on braneworlds was instigated by the 
proposals of Arkani-Hamed, Dimopoulos and Dvali 
(ADD) 1] and Randall and Sundrum (RS) (see for 
early work, and 4] for braneworld reviews). The ADD 
model has n flat, compact extra dimensions of size R. 
Due to the confinement of ordinary matter to a brane, R 
can be as large as ~ 0. 1mm (the scale down to which New- 
ton's law has been experimentally tested), and the model 
provides a possible resolution to the hierarchy problem if 
n > 2. More interesting from the general relativity view- 
point are the RS models proposed shortly afterwards. 
They allowed the bulk geometry to be curved, and en- 
dowed the branes with a tension. Their first model con- 
sisted of two branes of equal but opposite tension bound- 
ing a slice of Anti-de-Sitter space. This model also gives 
a possible resolution to the hierarchy problem, provided 
we live on the negative tension brane. In their second 
model, RS considered a single, positive tension brane in 
an infinite bulk. This model, which is loosely motivated 
by string theory |5( , has been one of the most popular to 
explore, and is the one we will be using. 

Even though the extra dimension is infinite and grav- 
ity is inherently five dimensional, RS showed that the 
Newtonian potential of a particle on the brane was in- 
deed the 4D 1/r potential. This result was backed up by 
more complete analyses which confirmed that the gravi- 
ton propagator did indeed have the correct tensor struc- 
ture, and that the effect of the extra dimension was to 
introduce a 1/r 3 correction to the gravitational poten- 
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An elegant description of non-perturbative gravity on the 
brane was provided by Shiromizu, Maeda and Sasaki 0. 
Using a Gauss-Codazzi approach, they projected the 5D 
Einstein equations onto the brane to obtain the effective 
4D field equations: 
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Here, A4 is a residual cosmological constant on the brane 
and represents the mismatch between the brane tension 
and the negative bulk cosmological constant. T^ v is the 
usual energy-momentum tensor of matter on the brane, 
and S^i^T 2 ) consists of squares of T M „ and thus is a lo- 
cal, high energy correction term. £^ v consists of the pro- 
jection of the bulk Weyl tensor onto the brane, and is 
non-local from the brane point of view. It is important 
to emphasise that since £ M „ is not given in terms of data 
on the brane, the system of equations J5J is not closed, 
in general. 

The generalisation of the FRW Universe that follows 
from Eq. (J2J has been well explored @. The term 
contributes a high energy correction term to the Fried- 
mann equation, which is relevant only in the very early 
Universe, and the Weyl term contributes a 'dark radi- 
ation' term. Hence, although all the implications may 
not have been calculated, braneworld cosmology for the 
pure RS scenario is pretty well understood. The situa- 
tion for braneworld black holes (BBHs) is somewhat more 
complicated however, and there is no longer a simple so- 
lution 0,^3- Black holes are fascinating objects, and 
provide a potential testing ground for general relativity. 
It is therefore important to investigate braneworld gener- 
alisations of the Schwarzschild solution, and the possible 
observational signatures that could result. 

The theory of gravitational lensing has been mostly 
developed in the weak field approximation, where it has 
been successful in explaining all observations ^l|. How- 
ever, one of the most spectacular consequences of the 
strong gravitational field surrounding a black hole is the 
large bending of light that can result for a light ray pass- 
ing through this region. The study of strong gravitational 
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lensing was resurrected recently by Virbhadra and El- 
lis [lj , who studied lensing by the galactic supermassive 
black hole, in an asymptotically flat background. Frit- 
tclli, Kling and Newman [13| found an exact lens equa- 
tion without reference to a background metric and com- 
pared their results with those of Virbhadra and Ellis. 
In Q Bozza et al. first defined a strong field limit and 
used it to investigate Schwarzschild black hole lensing an- 
alytically. This technique has been applied to Reissner- 
Nordstrom black holes [H and the GMGHS O charged 
black hole of heterotic string theory 01 1 an d was gener- 
alised to an arbitrary static, spherically symmetric metric 
by Bozza |l8j . In this paper we utilise this method to in- 
vestigate the gravitational lensing properties of a couple 
of candidate BBH metrics. 

Similar studies have been performed for a BBH with 
the induced geometry of the 5D Schwarzschild solution: 
9tt — 9rr = 1 — r h/ r2 - Both weak field lensing and 
strong field lensing |2(j for this geometry have been stud- 
ied. However, this metric is only appropriate for black 
holes with a horizon size smaller than the AdS length 
scale of the extra dimension: < I < 0.1mm. Hence, 
this metric is not appropriate for investigating the phe- 
nomenology of massive astrophysical black holes. 



II. BRANEWORLD BLACK HOLES 

The general static, spherically symmetric metric on the 
brane can be written as: 

ds 2 = g^dx^dx" = A 2 (r)dt 2 - B 2 (r)dr 2 - C 2 (r)d^. 

(3) 

Clearly, this is not in the simplest gauge, as we can still 
choose our radial coordinate, r, quite arbitrarily. The 
vacuum brane field equations following from Eq. @ (with 
A4 set equal to zero) are 
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The solution of these equations requires the input of 
from the full 5D solution. In the absence of such a solu- 
tion, an assumption about £ M „ or g^ v must be made in 
order to close the system of equations. 

Several special solutions, making various assumptions 
about g^, have been presented in the literature. The 
first attempt at a BBH solution was the so-called black 
string solution of Chamblin et al. |9( , which consists sim- 
ply of the 4D Schwarzschild solution 'stacked' into the 
extra dimension. Unfortunately, this solution has a sin- 
gular AdS horizon and is unstable to classical perturba- 
tions • The assumption A 2 = 1/B 2 leads to the tidal 
Reissner-Nordstrom solution of Dadhich et al. [231 : 
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Unlike the standard Reissner-Nordstrom solution, the 
'tidal charge' parameter Q can take both positive and 



negative values. Indeed, negative Q is the more natu- 
ral since intuitively we would expect the tidal charge to 
strengthen the gravitational field, as it arises from the 
source mass M on the brane (see (2^ for further dis- 
cussion). This metric has the correct 5D (~l/r 2 ) short 
distance behaviour and so could be a good approximation 
in the strong field regime for small black holes. Solutions 
have also been found which assume a given form for the 
time or radial part of the metric [23ll2"ij . Visser and Wilt- 
shire [2^ presented a more general method which gener- 
ated an exact solution for a given radial metric form. 

In all the above cases the radial gauge C — r was 
chosen (although [25| did comment on how to use their 
method when C(r) was not monotonic). However, there 
are good reasons to believe that the area A of the 2- 
spheres might not be monotonic. The second derivative 
of the area radius C (i.e. the radial function defined by 
y/A/4n) is given by 
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Hence for the area function to be guaranteed to be mono- 
tonic we must have G\~ G r r > 0, which is equivalent to 
the dominant energy condition. While this is generally 
satisfied in standard Einstein gravity, it need not be in 
the case of extra dimensions, and so it is important for 
the exploration of BBHs that we do not make the restric- 
tive ansatz C = r. (For a discussion of non-monotonic 
radial functions in the context of braneworld wormhole 
solutions, see 26].) 

An alternative to making guesses for the metric func- 
tions g^y is instead to make an assumption about the 
Weyl term £^ u . Although is a complete unknown 
from the brane point of view, the s ymm etry of the prob 
lem allows it to be decomposed as 
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where 



is a unit time vector and r M is a unit radial 



vector. Recently, we proposed a pragmatic approach to 
BBHs, in which an equation of state for the Weyl term 
is assumed [28| (see also p^t: 
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Of course, a priori there is no reason to suppose that the 
Weyl term should obey an equation of state. However, 
it is quite possible that it might have certain asymptotic 
equations of state which may be useful as near-horizon 
or long-range approximations to the (as yet unknown) 
exact solution. Using a dynamical systems approach, the 
system of equations (@J was solved, and the behaviour 
of the solutions classified according to the equation of 
state parameter 7. It was found that asymptotically flat 
solutions require an equation of state 7 < 0, and for 
I7I > 3 the BBH solutions have a singular horizon, and 
are allowed both with and without turning points in the 
area function. 



Using holography considerations, it was argued that 
we might expect equations of state with large 7 to be 
relevant near the horizon. Taking this reasoning to its 
extreme, we proposed as a 'working' metric for the near- 
horizon geometry the hi = (i.e. 7 = ±00) solution: 
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which has a turning point in the area function at r = rt , 
and the horizon at r — is singular (except for the spe- 
cial case rh = r t , which is just the standard Schwarzschild 
solution in isotropic coordinates). This metric has ap- 
peared in the area gauge as [2i| 
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1 In the context of the equation of state, JSJ has 7 



D 03 



where R = (r+r t ) 2 /r, GM — 2r t and GMe = rh—r t . For 
e > this gauge is valid in the whole horizon exterior, 
however for e < the turning point r t is outside the 
horizon Th and so the area gauge is inappropriate. 

Although the choice of metric © is somewhat arbi- 
trary, we believe that the horizon is likely to be singular 
and that a turning point in the area function is also likely, 
and this metric exhibits both these features. It has the 
added advantage of being analytic, and so seems a good 
choice for exploring these radical differences to the stan- 
dard Schwarzschild geometry. In this paper we investi- 
gate the gravitational lensing properties of the metrics 
(0 1 and (JSJ, to see how braneworld effects might mani- 
fest themselves in observations of black holes. 

It is important to emphasise that we envisage these 
only as possible near-horizon asymptotes of a more gen- 
eral metric, which has yet to be found. Neither metric 
satisfies the long distance 1/r 3 correction to the gravi- 
tational potential, and both would be constrained in the 
weak field by the PPN observations. 



III. GRAVITATIONAL LENSING 

The lensing setup is shown in Fig. ^ Light emitted by 
the source S is deflected by the lens L and reaches the 
observer O at an angle 8 to the optic axis OL, instead of 
p. The spacetime, described by the metric (jSJ centred on 
L, is asymptotically flat, and both observer and source 
are located in the flat region. By simple trigonometry, 
the lens equation can be written down: 

tan/3 = tan 6 - — — [tan0 + tan(a - 6)) . (11) 
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FIG. 1: Gravitational lensing diagram. 



From the null geodesic equations it is straightforward to 
show that the angular deflection of light as a function of 
radial distance from the lens is 
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By conservation of angular momentum, the impact pa- 
rameter u is given by 
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where the subscript indicates that the function is eval- 
uated at the closest approach distance ro. Hence, the 
deflection angle is given by 



a{r ) = I(r ) - it 
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Equations i|ll|) and 114f) are the basic equations of grav- 
itational lensing. In principle, the deflection angle a for 
a given metric can be calculated from Eq. I|14|) . This can 
then be plugged into the lens equation (|llfl . and the im- 
age position 9 for a given source position P can be found. 

The theory of gravitational lensing has been developed 
mostly in the weak field limit, where several simplify- 
ing assumptions can be made. The angles in Eq. 1|11|) 
are taken to be small, so that tan a; can be replaced by 
x for x = P,8, a, and the integrand in Eq. (|14|> is ex- 
panded to first order in the gravitational potential. For 
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FIG. 2: General behaviour of the deflection angle as a function 
of ro- As ro decreases, a increases, and each time it reaches a 
multiple of 2tt the photon completes a loop around the black 
hole. 



the Schwarzschild geometry, and setting f3 
to the well-known result: 
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0, this leads 
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where 9e is the Einstein radius. In this formulation, 
general relativity has been successful in explaining all 
observations (see for detailed reviews). However, it 
is important that gravitational lensing is not conceived of 
as a purely weak field phenomenon. Indeed, gravitational 
lensing in strong fields is one of the most promising tools 
for testing general relativity in its full, non-linear form. 



A. Strong field limit 

As the impact parameter u of a light ray decreases, the 
deflection angle a increases as shown in Fig. [21 At some 
point, the deflection angle exceeds 2ir and the photon 
performs a complete loop around the black hole before 
emerging. The images thus formed are termed 'relativis- 
tic images' and a theoretically infinite number of such 
images are formed on either side of the lens, correspond- 
ing to successive winding numbers around the black hole. 
The photon sphere is the radius r p at which a photon can 
unstably orbit the black hole, and is defined as the largest 
solution to the equation 
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As r approaches r pi with corresponding impact param- 
eter 
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the deflection angle diverges and for ro < r p the photon 
is captured by the black hole. 



Bozza [l8j has shown that this divergence is logarith- 
mic for all spherically symmetric black hole metrics of 
the form JSJ) . Hence the deflection angle can be expanded 
close to the divergence in the form 
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or in terms of the angular position of the image, 9 
u/D i, 
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-1 )+b + 0(u-u p ), (19) 



where the strong field limit (SFL) coefficients a and b 
depend on the metric functions evaluated at r p . This 
formula allows a simple, analytic description of the rela- 
tivistic images and their properties, rather than having 
to use the exact deflection angle calculated numerically 
from Eq. {11} ■ 

Equation (|19fl can be derived from Eq. i|14|) by splitting 
the integral into a divergent and non-divergent piece, and 
performing some expansions. Defining the new variable 
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the integral in equation l|14l) becomes 
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where 
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The function R(z, ro) is regular for all values of z and ro, 
but f{z, ro) diverges for z — * 0. Expanding the argument 
of the square root in /(z, ro) to second order in z: 



f(z, ro) ~ fo(z, r ) = (m(r )z + n(r )z 2 ) 
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(26) 



it is clear why the deflection angle diverges logarithmi- 
cally for ro = r p : with r p given by Eq. I|16f) . w(r p ) 
vanishes. Hence for ro = r pi fo oc 1/z and the integral 
(|21|l diverges logarithmically. 

Proceeding to split the integral H21(l into a divergent 
and a regular piece, and performing further expansions 
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(see ^3 f° r the detailed derivation 2 ), the SFL coefficients 
are obtained: 
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and 



n P = n(r p ) = { \ c ^ K A 'p ~ C 'p A p) ( 29 ) 
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[R(z, r p )f(z, r p ) - R(0, r p )f (z, r p )] dz, (30) 



is the integral I(r p ) with the divergence subtracted. 



FIG. 3: SFL coefficients for the U = metric 131L as func- 
tions of rh- The standard Schwarzschild case is given by 
r h = 1/4. 



IV. BRANEWORLD BLACK HOLE LENSING 

In this section we apply the method outlined in the 
previous section to calculate the deflection angle in the 
strong field limit for the candidate BBH metrics JSJ and 
|jSJ discussed in Section HT1 

A. {7 = metric 

The {7 = metric (jHJ has a number of key differences 
to the standard Schwarzschild geometry; the horizon is 
singular, and the area function can have a turning point 
that lies either inside or outside the horizon. Another 
difference, due to the fact that gtt ^ Qrr i ^ s that the 
ADM mass and gravitational mass (defined by gtt) are 
no longer the same. 

Normalising the distances to Art (which corresponds 
to a distance 2GM where M is the ADM mass - see 
Eq. I|10|)l. the metric functions are 
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(31) 



The radius of the photon sphere is given by: 

r P = | ( 1 + 4r h + ^l + Ar h + \Qrj \ . (32) 



2 Note: The expressions in llgl look slightly different to ours since 
we define the metric as g M „ = diag(A 2 , — B 2 , — C 2 ) as opposed 
to c; M „ = diag(A, —B, -C) in 



The SFL coefficients a, b and it p , calculated from 
Eqs. and 03, are shown in Fig. It 

can be seen that the biggest deviation from standard 
Schwarzschild lensing is for the minimum impact param- 
eter u p . This is because as the horizon is shifted in- 
wards/outwards relative to the Schwarzschild case, the 
photon sphere is pulled/pushed along with it. 

We can check the accuracy of the SFL approximation 
by comparing the exact deflection angle a GX act calculated 
from Eq. (H) with the SFL a S FL from Eq. jf|j. The out- 
ermost relativistic image appears where a ~ 27r, which 
occurs for an impact parameter u\ — u p + x, where 
x ~ 0.003 depends on r/j. The discrepancy between 
acxact(wi) and asFL(wi) is less than 0.13% for all val- 
ues of rh we consider. Hence, the SFL of the deflection 
angle is very accurate and can be reliably used to ob- 
tain accurate results for the properties of the relativistic 
images. 



B. Tidal Reissner-Nordstrom metric 

The tidal RN metric JSJ has the same properties as the 
standard RN geometry for q > 0: there are two horizons, 
both of which lie within the Schwarzschild horizon, and 
the singularity at r = is timelike. However, we can 
now have q < in which case there is just one horizon, 
lying outside Schwarzschild, and the central singularity 
is spacelike, as in the Schwarzschild case. Normalising 
the distances to 2GM, with q = Q/(2GM) 2 , the metric 
© becomes 



B 2 (r) = (l-]; + ^y\ 03) 
C 2 (r) = r 2 . 
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FIG. 4: SFL coefficients for the Tidal Reissner-Nordstrom 
metric 1331 . as functions of q. The standard Schwarzschild 
case given by q = 0. 



The radius of the photon sphere is given by: 



(34) 



The SFL coefficients are shown in Fig . 0] These re- 
sults reproduce those of Eiroa et al. 15] for q > (see 
also Hg), but we have extended the results to negative 
q. We emphasise that there is no electric charge present 
for the tidal RN solution - q is a tidal charge parameter 
arising from the bulk Weyl tensor. 

Again, comparing a CX act with osfl for an impact pa- 
rameter corresponding to the outermost image, it is 
found that the discrepancy is less than 0.5% for the val- 
ues of q considered here. 



V. OBSERVABLES 



In Section IlII Al it was shown how to calculate the de- 
flection angle in the strong field limit and in Section HVI 
this method was applied to the candidate BBH metrics. 
In this section we put the SFL of the deflection angle 
into the lens equation to obtain analytic formulae for the 
properties of the relativistic images in terms of the SFL 
coefficients a, b and u p . 

As expected, the relativistic images formed by light 
rays winding around the black hole are greatly de- 
magnified compared to the standard weak field images, 
and are most prominent when the source, lens and ob- 
server are highly aligned |l2j| . Hence, we restrict our 
attention to the case where and 9 are small (see [3(j 
for the general case where this assumption is relaxed). 
Although we can not assume a is small, if a light ray 
is going to reach the observer after winding around the 
black hole, a must be very close to a multiple of 2tt. 
Writing a = 2nir + Aa n , n £ Z, the lens equation 
becomes 
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Firstly, we have to find the values 6„ such that a (6*° ) 
2nir. With a given by Eq. (|19|l we find 



9 = J^P/j 
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where 
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Thus the position of the n th relativistic image can be 
approximated by [is| 
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where the correction to 6„ is much smaller than 0„ . Ap- 
proximating the position of the images by 0°, the mag- 
relativistic image is given by 
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Equations (|38|l and 1)39(1 relate the position and magnifi- 
cation of the relativistic images to the SFL coefficients. 
We now focus on the simplest situation, where only the 
outermost image 9\ is resolved as a single image, with 
the remaining images packed together at 9^ = u p /D i. 
Therefore we define the observables 



s = 
I = 



Mi 



(40) 
(41) 



which are respectively the separation between the outer- 
most image and the others, and the flux ratio between 
the outermost image and all the others. It is found that 
these simplify to |18|: 



s 
f 



J2-K/S. 



(b-2ir)/a 



(42) 
(43) 



These equations are easily inverted to give a, b and so 
if an observation were able to measure s, / and 9^ the 
SFL coefficients could be determined and the nature of 
the lensing black hole identified. 



An Example: The galactic supermassive black 
hole 



It is believed that the centre of our galaxy harbours a 
black hole of mass M = 2.8 x 10 6 M© [3l|. Taking D ol = 
8.5 kpc, Virbhadra and Ellis 01 studied the lensing of a 
background source by this black hole and found that the 
relativistic images are formed at about 17 /i arc sec. from 
the optic axis. 

In Table [I] we estimate the observables 9^, s, f defined 
in the previous section for the U = and tidal RN BBH 
metrics, as well as the standard Schwarzschild metric. It 
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Schwarzschild 


17 = 


Tidal RN 


r h = 0.1 


r h = 0.2 


r h = 0.3 


r h = 0.4 


q = -0.2 


q = -0.1 


q = 0.1 


q = 0.2 


9oo (ft arc sec.) 
s (n arc sec.) 
f m (mags.) 


16.87 
0.0211 
6.82 


13.24 
0.0303 
6.13 


15.65 
0.0235 
6.61 


18.11 
0.0192 
7.01 


20.62 
0.0164 
7.32 


18.85 
0.0150 
7.26 


17.92 
0.0173 
7.08 


15.64 
0.0286 
6.44 


14.07 
0.0502 
5.70 



TABLE I: Estimates for the lensing observables for the central black hole of our galaxy. Ooa and s are defined in Section [V] 
and fm ~ 2.5 log / is / converted to magnitudes. 



is clear that the easiest observable to resolve is 9^, since 
a microarcsecond resolution is in principle attainable by 
VLBI projects such as MAXIM [H and ARISE [if. 
However, the disturbances inherent in such observations 
would make the identification of the faint relativistic im- 
ages very difficult, as discussed in [T^|. 

If a measurement of 9^ was made, it would be imme- 
diately capable of distinguishing between Schwarzschild 
and other types of geometry. However, to determine all 
the SFL coefficients and thus unambiguously identify the 
nature of the lensing black hole, it is necessary to also 
measure s and /. This would require the resolution of 
two extremely faint images separated by ~ 0.02 \i arc sec. 
Such an observation in a realistic astrophysical environ- 
ment is certainly not feasible in the near future, although 
if such an observation were ever possible, it would provide 
an excellent test of gravity in a strong field. 

VI. CONCLUSIONS 

Gravitational lensing in strong fields provides a poten- 
tially powerful tool for testing general relativity, and the 
strong field limit provides a useful framework for com- 
paring lensing by different black hole metrics. Of the 
possible alternatives to standard GR, braneworld grav- 
ity is a very interesting model to explore given the current 
interest in theories with extra dimensions. 

In this paper we have investigated strong field lensing 
by potential near-horizon BBH metrics. Although the 
correct BBH metric is unknown and much theoretical 
work remains to be done, this study is a useful first step to 
explore the possible effects that braneworlds could have 
on the spacetime surrounding a black hole. 

Table [I] clearly shows that BBHs could have signifi- 



cantly different observational signatures than the stan- 
dard Schwarzschild black hole. Although the resolutions 
required for these observations are beyond reach of cur- 
rent observational facilities, this encourages the inves- 
tigation of more realistically observable situations. An 
interesting possibility in this direction is the study of the 
accretion discs surrounding black holes. 

The observed disc emission depends on several factors 
that could get modified by braneworld effects. A key 
factor is the radius of the innermost stable circular orbit, 
since emitting material at this radius sets the maximum 
temperature for the disc emission |34| . Just as the ra- 
dius of the photon sphere is shifted inwards or outwards 
relative to the Schwarzschild case for a BBH, so too is 
the radius of the innermost stable orbit for matter pgj . 
In addition, the observed disc emission, and in particu- 
lar the iron fluorescence line profile, is affected by rela- 
tivistic effects [3j| (doppler shift, gravitational redshift), 
which would be modified if the metric in the emitting re- 
gion was not that of standard general relativity, but was 
a modified braneworld metric. Also, the light rays are 
gravitationally lensed by the central black hole as they 
escape the disc, and we have shown here that such lensing 
can be different from the Schwarzschild case for BBHs. 
In light of the results found here, it is not unreasonable 
to anticipate that these effects could result in distinctive 
observational signatures for accretion discs. This work is 
in progress. 
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